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Observation of spin-selective tunneling in SiGe nanocrystals 

G. Katsaros, 1 ^ V. N. Golovach, 1 P. Spathis, 1 N. Ares, 1 M. Stoffel, 2 
F. Fournel, 3 0. G. Schmidt, 2 L. I. Glazman, 4 and S. De Franceschi 1 

1 SPSMS, CEA-INAC/UJF-Grenoble 1, 17 Rue des Martyrs, 38054 Grenoble Cedex 9, France 

2 Institute for Integrative Nanosciences, IFW Dresden, Helmholtzstr. 20, D-01069 Dresden, Germany 

3 CEA, LETI, MINATEC, F38054 Grenoble, France 

4 Department of Physics, Yale University, New Haven, Connecticut 06520, USA 

(Dated: July 21, 2011) 

Spin-selective tunneling of holes in SiGe nanocrystals contacted by normal-metal leads is reported. 
The spin selectivity arises from an interplay of the orbital effect of the magnetic field with the strong 
spin-orbit interaction present in the valence band of the semiconductor. We demonstrate both 
experimentally and theoretically that spin-selective tunneling in semiconductor nanostructures can 
be achieved without the use of ferromagnetic contacts. The reported effect, which relies on mixing 
the light and heavy holes, should be observable in a broad class of quantum-dot systems formed in 
semiconductors with a degenerate valence band. 



The spin-orbit interaction (SOI) has become of cen- 
tral interest in the past years^, because it enables an all- 
electrical manipulation of the spin. In the field of spin 
qubits, one of us 2 - suggested the electrical control of local- 
ized spins by means of the electric-dipole spin resonance, 
and this scheme has been successfully used for spin ro- 
tations of electrons in quantum dots (QDs) 3,4 . Already 
much earlier, Datta and Das 5 proposed a semiconductor 
transistor that would operate through a gate-controlled 
spin precession, mediated by the SOI. In this type of spin 
transistor, spin-polarized electrons are injected into the 
semiconductor from a ferromagnetic (FM) contact. The 
realization of an efficient spin injection has proven to be 
a difficult task 6,7 . Only recently, high spin-injection ef- 
ficiencies were reported for FM contacts to semiconduc- 
tors^—. In nanostructures, however, experimental ev- 
idence of spin injection is not as strong and clear— ~—. 
Here we show that the SOI in the valence band, quan- 
tified by the spin-orbital splitting Aso, provides an al- 
ternative way to obtain spin-selective tunneling without 
requiring FM electrodes. 

At cryogenic temperatures, transport through QDs is 
dominated by the Coulomb blockade (CB) effect. In the 
CB regime, single-hole transport is suppressed and elec- 
trical conduction is due to second-order cotunneling (CT) 
processes — . We consider here the case of a QD with an 
odd number of holes and a spin-doublet ground state. 
A magnetic field, B, lifts the spin degeneracy by the 
Zeeman energy Ez = gfisB, where g and fis are the 
hole g-factor and Bohr magneton, respectively. Once the 
bias voltage across the QD exceeds the Zeeman energy, 
\eV\ > Ez, the inelastic CT processes can flip the QD 
spin, leaving the QD in the excited spin state; hereinafter 
e is the elementary charge (e > 0). The onset of spin-flip 
inelastic CT manifests itself as a step in the differential 
conductance, G = dl/dV, at eV = ±Ez 17 ■ 

Our measurements reveal a pronounced asymmetry in 
the step height of G with respect to the polarity of V, as 
recently predicted by Paaske et a/. 1? , in a model with a 
rather generic form of the SOI interaction. The asymme- 
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FIG. 1. (color online) Spin-selective tunneling in (a) a QD 
coupled to FM leads and (b) a QD with SOI coupled to 
non-magnetic leads. The solid (dashed) arrows indicate the 
tunneling processes involved in the inelastic CT for forward 
(reverse) biasing, with solid (dashed) arrows representing 
stronger (weaker) tunnel rates. In both (a) and (b), the tun- 
nel rate, V = ttv \t\ , differs for each Zeeman sublevel of the 
QD. In setup (a), it is the density of states v that brings 
about the spin selectivity of the tunneling. In setup (b), the 
spin selectivity is caused by the tunneling amplitude t, which 
is sensitive to the spinor wave functions at the point of tun- 
neling. In the valence band, for energies < Aso, the B-field 
efficiently makes T spin-dependent by affecting the mixing be- 
tween heavy and light holes. Since the inelastic CT current 
is proportional to T^rJ for the forward bias and to T^rJ for 
the reverse bias, an asymmetric G(V) is expected whenever 
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try is found to depend on the magnitude and direction of 
B. Our results are consistent with an explanation based 
on the Luttinger Hamiltonian for the valence band of the 
semiconductor. The spin selectivity of tunneling arises 
from an interplay of the complex structure of the valence 
band with the orbital effect of the magnetic field. At 
B = 0, the time-reversal symmetry ensures that the two 
states forming the Kramers doublet in the QD are in- 
distinguishable and the spin selectivity of tunneling van- 
ishes. 

It is interesting to note that the transport characteris- 
tics of a QD with SOI coupled to normal leads are similar 
to those of a QD without SOI coupled to FM leads. We 
illustrate this similarity in Fig. [TJ where we consider the 




FIG. 2. (color online) (a) Schematic of a QD device fab- 
ricated from a SiGe self-assembled nanocrystal grown on a 
silicon-on-insulator substrate having a heavily doped handle 
wafer which is used as a back gatei^. (b) Qualitative band 
diagram of a Ge-rich SiGe quantum well illustrating the effect 
of quantum confinement along the growth (z) direction: HH 



and LH branches are split at k x 
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and anti-cross at 



finite k x or k y . The red dots indicate that many other HH 
subbands exist before the first LH subband is encountered, 
(c) G(V) for different perpendicular B-fields from 0.5 to 3 T. 
Inset: G(Vg, V) for a 75-mT perpendicular field (Vg spans a 
range of 850 mV and V ranges from -3.5 to 3.5 mV). (d) G(V) 
for different parallel fields from 1 to 7 T. The Zeeman split- 
ting of the Kondo peak is asymmetric in (c) and symmetric 
in (d). (e) Angular dependence of the split Kondo peak for a 
fixed Vg and B — 3T. Superimposed G(V) traces for 8 — 0, 
30, 60, and 90 degrees. 



simplest case, in which the Zeeman interaction and the 
two spin-dependent tunnel contacts have collinear quan- 
tization directions. 

We have studied the low-temperature magneto- 
transport properties of individual SiGe self-assembled 
QDs with a base diameter d « 80 nm and a height 
w ~ 20 nm. The hole motion is strongly quantized along 
the growth direction [001]. A schematic of a typical QD 
contacted with Al electrodes is shown in Fig. [2ja). For 
such QDs, the hole wave function is generally composed 
of both heavy holes (HHs) and light holes (LHs). Due 
to the confinement and compressive strain, the degener- 
acy between the HH and LH branches, present in bulk 
at the T-point, is lifted. The HHs become energetically 
favorable. In Fig. [2jb), we illustrate the interaction be- 
tween a HH and a LH branch in the two-dimensional 
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FIG. 3. Asymmetry parameter A as a function of (a) 9 and 
(b) perpendicular B. We note that, by subtracting the elastic 
CT contribution from G(V), A becomes larger than 0.5. 



(2D) case. The split-off band is far away in energy due 
to a large Aso ■ HHs and LHs are states of angular mo- 
mentum | with projections ±| and ±A, respectively, see 
Appendix IB] We remark that, contrary to the HH states, 
the LH states cannot be factorized into a product of spin 
and orbital components. 

The stability diagram, G(Vg,V), of a QD device is 
shown in the inset of Fig. H{c). The diamond-shape re- 
gion delimited by dashed lines highlights the CB regime 
for an odd number of confined holes. While G is gener- 
ally suppressed within this CB diamond, a G resonance 
can be identified at V — 0, providing a clear signature of 
a spin-1/2 Kondo effect*^. At finite B, this resonance is 
split by the Zeeman effect as shown in Figs. [2] (c) and (d) 
for perpendicular and parallel B, respectively (all G(V) 
traces were taken at the same Vg). 

For perpendicular B [Fig. Efc)], the splitting of the 
Kondo peak is clearly asymmetric with respect to a sign 
change in V. The asymmetry in G arises at the onset of 
spin-flip inelastic CT (i.e. for \eV\ > Ez). For parallel B, 
however, the asymmetry is practically absent [Fig. |2fd)]. 
To further investigate this anisotropy, a sequence of G(V) 
traces was taken while rotating a 3 T field in a plane per- 
pendicular to the substrate. The resulting data, G(9, V), 
are shown in Fig. HJe), with 9 being the angle between 
the field and the substrate plane. Along with a variation 
in the Zeeman splitting of the Kondo peak, caused by 
the ^-dependent hole g-factor— , the asymmetry becomes 
progressively more pronounced when going from 9 = 
(or 180°) towards Q = 90°. 

The asymmetry observed in G(V) can be quantified 
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where G± = G(±E z /e). The detailed 



G-+G+ 

A{9) dependence, extracted from Fig. ^e), is shown in 
Fig. ^a). A w for 9 = (or 180°) and it increases 
monotonically up to 0.2 for 9 approaching 90°. The 
same qualitative behavior was observed in a second de- 
vice, which did not display Kondo effect, see Appendix [Al 
The asymmetry A reaches 0.4 at 3 T for that device. We 
remark that, although the first device shows larger con- 
ductance due to the Kondo effect, the asymmetry A in 
both devices is a consequence of spin-dependent tunnel 
rates. 

In order to explain the microscopic origin of the mea- 
sured effect, we represent the Luttinger HamiltonianSi as 



a block matrix in the basis of HH (h) and LH (I) states, 

Hhh U hi 



H = 



Hm Hn 



(1) 



In the blocks Hhh and Hu , we discard all terms that van- 
ish in the 2D limit (w/d — > 0), whereas in the blocks Hhi 
and Hih, we keep only the leading-order in w/d terms. 
A systematic expansion around the 2D limit is outlined 
in the Appendix [B] The blocks Hhh and Hu assume a 
familiar form 
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+ ^bo- ■ g h/ i ■ B + U(x, y) + V h /i {z), (2) 

where the axes x, y, and z point along the main crystal- 
lographic directions, with z = [001] being the direction of 
the strongest quantization. After the expansion around 
the 2D limit, the kinetic momentum operators k x and 
k y contain only the component B z , whereas k z is inde- 
pendent of B. In Eq. @ and below, 71, 72, 73, k, and 
q denote the Luttinger parameters-- and m denotes the 
bare electron mass. The Pauli matrices <x = (a x , a y , a z ) 
represent the remaining pseudo-spin degree of freedom in 
each block. We choose the following pseudo-spin basis 30 : 



|t) h = |3/2,-3/2), \{) h 
|t>, = |3/2, +1/2) , \i) r - 



: |3/2,+3/2), 
13/2,-1/2). (3) 



In the (x, y, z)-frame, the tensors of the g-factor are di- 
agonal: gh = diag (0, 0, — 6k) and gi = diag (4k,4k,2k), 
where we neglected, for simplicity, the terms proportional 
to the smallest Luttinger parameter q. The minus sign 
in (gh)zz is due to our basis choice in Eq. ^. In Eq. @, 
we included an in-plane confining potential U(x,y). The 
motion along z is confined to an infinitely-deep square 
well, with different offsets, Vh and Vi > Vu, due to strain. 
The blocks Hhi and Hih are given by 



+ v 373 

H h i = (Hih) = i (k x a y + k y a x ) k z 



(4) 



These blocks intermix HHs and LHs, such that the wave 
function of the hole in a given QD state assumes the 
general form W = a^h + P^i- In terms of the true-spin 
states, such a wave function consists of a superposition of 
the spin-up (j) and spin-down (1) states entangled with 
the orbital degrees of freedom: 



%W= $ iWt+XiW; 
%(r)=xa(r)t+*2(r)+, 



(5) 



where ff and J| denote the components of the Kramers 
doublet in the QD. Focusing on the first HH subband, 
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FIG. 4. (color online) (a) Tunneling polarization M z as a 
function of the space coordinate p for two Fock-Darwin states 
(n, to) as indicated and for uj c = O.lwo- (b) _B-field depen- 
dence of the difference AM Z = M z (pl) — M z (pn), for the 
values of pL and pn indicated in (a) . (c) AM Z as a function 
of 9 for a fixed value of \B\. The value of \B\ corresponds to 
loc — wo at 9 = 90° . 
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and similar expressions for $2(0 and X2( r ), obtained 
from Eq. ^ by replacing U- — > U +1 fc_ — > k +1 
and iph{x,y) — > ^hfoiV)- ^ n our notation, U± = 
T-ys (A ± iY) and k± = T-js (k x ± ik y ). The Bloch am- 
plitudes A, Y, and Z describe the valence band in the 
absence of SOI. In blocks Hhh and Hu, the motion along 
z separates; we denote the corresponding eigenenergies 
and eigenfunctions by E n and f n (z), respectively. 
The tunneling amplitudes t\ s are found as 



u=X,Y,Z 



T u (u,<r\V a (n)), 



(7) 



where T u is the coupling strength between Bloch ampli- 
tude u and the lead, and (u, a\^ s (ri)) are the projec- 
tions of the QD eigenstates $" s (r), see Eq. ([5]), onto the 
product state of Bloch amplitude u and spinor \o). The 
tunneling amplitudes in Eq. ([7]) depend on the point of 
tunneling, t\ = tl, r^, the component of the true spin in 
the lead, a =t i, and the component of the Kramers dou- 
blet on the dot, s =f|-, JJ.. We remark that T X ,T Y , and T z 
appear in Eq. J7J as phcnomcnological parameters. They 
depend on the details of the metal-semiconductor inter- 
face and cannot be determined within the k • p-theory 
used here. We find 
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where $j(r) and Xi( r ) are obtained from Eq. ([6]) by re- 
placing Z -> T z and U± -> T^j {T x ±iT Y ). 

The spin selectivity of the tunneling is best seen in the 
matrix of the tunnel rates, T ss / = 7r ^ CT t* as v a t as i. With 



v\ = v\, (case of non-FM leads) we find, up to a common 
factor, 
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At B = 0, time-reversal symmetry requires that 
*a(r) = [*i(r)]* and xa(r) = - ttt.(r)]* , 



(9) 



(10) 



leading to T^ = T^ and r^ = T^ = in Eq. ([9]). 

At B 7^ 0, however, the orbital effect of the i?-ficld 
modifies the functions $j(r) and Xi( r )j such that the re- 
lations in Eq. (|10[) are no longer satisfied. In general, the 
matrix T ss > has nonzero off-diagonal elements. Since it 
is a hermitian matrix, there exists a direction in space, 
M, such that an SU(2) rotation of the Kramers dou- 
blet by an angle Z.zM makes the rate matrix diagonal, 
r = diag(r^, r^), with T^ > r^. To quantify the spin 
selectivity of the tunneling, we define 



\M\ = 



r t ~ r 4J 

r^ + r,. 



(ii) 



In respect to transport, M is analogous to the polar- 
ization vector of the FM. Indeed, the maximum of spin 
selectivity in tunneling from a FM is achieved when the 
FM is a half-metal, e.g., v^ ^ and v^ = 0. This extreme 
case corresponds to M — 1 and can be approached in our 
case by increasing B z . 

In order to illustrate the origin of the spin selectiv- 
ity, we focus on the special case: Tx — Ty — and 
Tz 7^ and refer to this tunneling model as the Z- 
model. In the Z-model, vector M is parallel to the z- 
axis. Tunneling to the hole states is possible only due 
to the admixture of the LH subbands. Furthermore, 
in this model, the spin selectivity is determined by the 
fact that_Xi(r) oc k^ip h (x,y) and X2(r) oc k + ipl(x,y), 
whereas &i(r) = 0. Using this information in Eqs. (|9]) 
and (TTTj) . we specify ?/>/i(x, y) to the Fock-Darwin states^. 
Therefore, we assume that U(x,y) in Eq. @ is given by 
U(p) = to*WqP 2 /2, where m* is the effective mass for in- 
plane motion, ujq is the oscillator frequency of the har- 



monic potential, and p = x 



y 



For the first two 



states (n = and m = 0, — 1), we obtain 



with f(p) = 2h/(m*ujp 2 ) - 1. Now M z depends both 
on B z and on p, see Fig. Eta). The spin polarization of 
two contacts positioned arbitrarily on a QD may differ 
significantly from each other, see, e.g., points p^ and 
Pr in Fig. SJa). The asymmetry in the inelastic CT is 
related to AM Z = M z {p L )-M z {p R ) ^ 0. AM Z increases 
with B z [Fig. @] (b)], displaying at the same time strong 
dependence on the S-field direction [Fig. @] (c)], in good 
qualitative agreement with the results in Fig. [3j 

The described joint effect of SOI and Zeeman splitting 
explains our experimental findings. In addition, it opens 
the door to an original scheme for measuring Rabi spin 
oscillations in hole confinement QDs. Let us consider 
a spin-1/2 QD in the CB regime under a perpendicular 
B of the order of a few T. In such a case, a transport 
characteristic of the type shown in Fig. 2(c) is to be 
expected. For V — 0, no current flows through the QD. 
Yet we suggest that a finite current could be generated 
by a resonant rf field (at frequency / = Ez/h) capable 
of inducing coherent oscillations between the Zeeman- 
split states of the QD. In fact, as the excited ff state 
becomes populated it can decay to the ground J| state by 
an inelastic CT process — a hole tunnels out of the QD 
from the ff state being replaced by another hole tunneling 
into the -IJ. state. Because ft and J| states have tunnel 
couplings with opposite asymmetries, in a configuration 
such as the one depicted in Fig. 1(b) the most favorable 
CT relaxation process would involve the transfer of a hole 
from the right to the left contact. Hence a net dc current 
could be driven by a continuous resonant irradiation. In 
addition, combining rf bursts with syncronized Vq pulses 
may enable the coherent control of the QD pseudo-spin 
states. In this scheme, well-defined pseudo-spin rotations 
would be performed in the deep CB regime (i.e. during a 
negative Vq pulse), whereas pseudo-spin read-out would 
take place in the CT regime. 
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age used in Fig. [2(a). The work was supported by the 
Agence Nationalc dc la Recherche (through the ACCESS 
and COHESION projects), US DOE Contract No. DE- 
FG02-08ER46482 (Yale), and the Nanosciences Founda- 
tion at Grenoble, France. G.K. acknowledges support 
from the Deutsche Forschungsgemeinschaft. 
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where lu = \/i>Jq + u^/A, and u> c — eB z /m*c. For these 
states, M z depends on B z but not on p, see Fig.[4|a). For 
B z ^ 0, the contacts will exhibit spin-dependent tunnel 
rates with the same polarization value M z regardless of 
the point-tunneling position. In such a case, no asymme- 
try in the inelastic CT is expected. 

The situation changes starting from n = and m = 1, 
where 



M z = 



uj u c 1 

— /(P) + T- TFT 

w c 4a; / (p) 



(13) 



Appendix A: Second device 

Similar asymmetries as the ones described in the main 
text were observed also for a second device, less strongly 
coupled to the metallic leads. Figure[5](a) is a plot of the 
differential conductance G versus the angle of magnetic 
field and the bias voltage V for a fixed value of 3 T. 
Differently to the first device the minimum g-factor is not 
observed for a parallel magnetic field but there is a shift 
by about 15-20 degrees. Some characteristic traces of G 
vs V taken at 20 (blue), 45 (red) and 90 (black) degrees 
are shown in Fig. [5] (b). Interestingly, the position of the 
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FIG. 5. (a) Evolution of the differential conductance dl/dV vs 
the angle of magnetic field and V. The amplitude of the mag- 
netic field is fixed to 3T. (b) Characteristic traces of dl /dV vs 
V for 20 (blue), 45 (red)and 90 (black) degrees, respectively. 
The traces have been shifted by Q.06x2e 2 /h for clarity, (c) 
Plot showing the evolution of the asymmetry A vs the angle 
of the magnetic field, (d) Plot showing the evolution of A vs 
the value of the perpendicular field. 



minimum g-factor coincides with the case for which no 
asymmetry is appearing (15-20 degrees). We believe that 
both observations are because the 2D plane of the wave 
function is not parallel to the substrate. 

The asymmetry A follows the same trends as were ob- 
served in the first device. It is almost zero for 15-20 
degrees (the position of the minimum g») and it obtains 
values of about 0.35-0.4 at large out of plane angles of the 
magnetic field [Fig. [5] (c)]. The difference in the value of 
A for positive and negative magnetic fields is attributed 
again to the different angle the 2D hole wavefunction 
plane forms with the magnetic field. Finally, Fig. [3] (d), 
verifies that the asymmetry increases with B z . 

We remark that in the absence of misalignment the 
asymmetry A obeys the relation 



A(B Z ) = A(-B z 



(Al) 



which holds within the experimental accuracy for the 
device described in the main text. Equation (|A1[) can 
be understood as follows. On the one hand, the spin- 
selective part of r ss / is proportional to the orbital B z 
and therefore it changes sign when flipping the direction 
of the magnetic field. On the other hand, the Zeeman en- 
ergy is also changing sign when flipping the direction of 
the magnetic field, swapping thus the roles of the ground 
and the excited state. As a result, the measured cotun- 
neling asymmetry A does not change sign when changing 



Appendix B: Expansion around the 2D limit 



Our starting point is the Luttingcr Hamiltonian 24 , 
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where m is the mass of the electron in vacuum, 71, 72, 73, 
k, and q are the Luttinger parameters, k — (k x ,k y ,k z ) is 
the momentum of the band electron, 



k 



' ih ir + - A ( r ) 

or c 



(B2) 



with A(r) being the vector potential due to the magnetic 
field, e the elementary charge (e > 0), and c the speed of 
light. Further, {. . .} denotes the symmetrized product, 
e.g. 



\k x kyf — — (k x ky + kyk x ) . 



(B3) 



and J = (J x , J y , J z ) are 4x4 matrices representing the 
spin J = 3/2 in a basis of choice. We choose the basis 2 ^ 



where 
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V2 



(X + iY) t, 



V2 



(X - iY) I, 



-= [- (X + iY) I +2Z t] , 



(B5) 



Here, the Bloch amplitudes X, Y, and Z are chosen to 
be real. They belong to the representation T25< of the 
group Oh and transform under the group operations as 
X = yz, Y = xz, and Z = xy^— . The basis functions 
X, Y, and Z form a subspace that is equivalent to the 
space of the angular momentum / = l 24 . The states 
in Eq. (JB5I) originate from the addition 2 ^ J = I + S, 
where S is the electron spin (S = 1/2) in the usual basis 
{t,|}- The subspace of J = 3/2 is shown in Eq. (1B5J) . 
whereas the subspace of J = 1/2 is neglected, because 
it corresponds to the split-off band, i.e. the band that is 
shifted, due to the spin-orbit interaction, by an amount 
Ago below the top of the valence band. The Luttinger 
Hamiltonian describes the very top of the valance band, 
at energies E <§; Ago- I n the basis given by Eqs. (|B4|) 



and (|B5[) . the matrices of J — 3/2 read 25 
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In Eq. (|B1|) . the axes x, y, and z are fixed along the 
main crystallographic directions of the cubic crystal. We 
choose the axis z to point along the growth direction 
of the nanocrystal, making it the axis of the strongest 
size quantization. The idea of expanding the Luttinger 
Hamiltonian around the 2D limit consists in regarding 
quantities like z and d/dz as proportional to w and 1/w, 
respectively. Here, w is the width of the 2D layer (i.e. 
height of nanocrystal), which is considered to be much 
smaller than the nanocrystal diameter d. On the other 
hand, quantities like x and d/dx are regarded as pro- 
portional to d and 1/d, respectively^. Before expanding 
in powers of w/d -C 1, it is convenient to represent the 
Luttinger Hamiltonian in a block form. 

We use two projection operators, ph and pi, which 
project on the subspaces of the heavy (h) and light (I) 
holes. In terms of J 2 , they are written as 



Ph 



Pi 



1 
2 V ' " 4 

1/9 -J? 



J: 



2 \4 



(B9) 



Ph and pi resolve the unity, Ph+Pi = 1, and have the usual 
properties of projection operators: p\ = Ph, pf = Pi, and 
PhPi = PiPh = 0. The Luttinger Hamiltonian in Eq. (|B1[) 
can then be written as follows 

H = (Ph+Pi)H(p h + Pl ) = H hh \h) (h\+H hl \h) (l\ 
+ H lh \l)(h\+H n \l)(l\, (BIO) 

which makes up a matrix in the (h, Z)-space, 



H = 



Hhh Hhi 
Hih Hu 



(Bll) 



Each element H^ can be represented as a 2 x 2 matrix in 
the space of the pseudo-spin. The blocks on the diagonal 



read 



Hhh = 



Hu 



(fc 2 +fe 2 ) + 



7i +72 
2m 

+ o W?c ■ 9h ■ B, 



71-272,2 
2m 



7i -72 



(kl + k 2 y ) + 



7i + 2 7: 



2 ,.2 



2m 
-oMbt ■ gi- B, 



2m 



(B12) 



where the g-factors gh and gi are tensors. In the frame 
(x,y,z), they are diagonal: 



9h 



3q 
3q 



-6k- ^-q 



(B13) 



and 



91 



4k + 10g 
4«+10g 








2«.+ i<7 



(B14) 



The off-blocks are related to each other by hermiticity, 



Hu = {Hi h ) 



(B15) 



For Hih, we have 
Hm 



i {{k x k z } G y + {k y k z } a x ) 



2 m (^ - k y) + l m {k x k y } (Jz 
+V3/zb f K + -q J ((T X B X - OyBy) . (B16) 

All we have done so far was to rewrite the Luttinger 
Hamiltonian in a block form. Next we proceed with the 
expansion in powers of w/d <C 1 as explained above. 
We allow for gauges of the form 



A x = zB y - ryB z , 



Ay — zB x 
A z =0, 



(1-t)xB z , 



(B17) 



where r is a real number expressing the remaining gauge 
freedom in two dimensions. After taking the 2D limit, 
we will be able to use a reduced (2D) vector potential, 
a = (a x ,a y ), which is given by 



a x = -ryB z , 
Oy = (1 - t)xB z 



(B18) 



Having in mind such a transition, we pull out the z- 
dependence from k x and k y , 



k x ^ k x -\ z, 

c 

b b - eBx 

it'll * ft/11 Al ■ 



(B19) 



Here, on the right-hand side, k x and k y do not depend 
on z anymore, because they are given in terms of the 2D 
vector potential a(x,y) as 



fox 



-ih 
-ih 



d_ 

dx 
d 



dy 



c 

e 
+ -a v . 



(B20) 



The next step is to substitute Eq. (JB19I) into the Lut- 
tinger Hamiltonian and to group the terms according to 
their order of w/d. 

The substitution of Eq. (IB19I) in the blocks Hhh and 
Hu produces linear in k x and k y terms which are not mul- 
tiplied by any Pauli matrix. Such terms can be gauged 
away after integration over z, since they correspond to 
a constant shift in a x and a y . They also admix higher 
heavy-hole subbands and slightly renormalizc the inplane 
mass, but this admixture, as well as the mass renormal- 
ization, vanishes in the limit w — > 0, because the corre- 
sponding perturbation is proportional to z. Therefore, 
we dispense with the new terms generated in Hhh and 
H u . 

We note that, for the blocks Hhh and Hu, the transi- 
tion to 2D is identical to what is usually done for electrons 
in the conduction band. 

For the blocks Hih and Hhi, we make the substitution 
in Eq. (JB19I) and obtain lots of terms. The origin of each 
term can be traced back through the following interme- 
diate step: 

eB 

\k x kzf ~ ^ k x k z H \zk z \ , 

c 

eB x 

\kykzj y kyk z \zk z \ , 



hi 



2e 
k x k -\ [k x Hy + kyli x ) z 



.y , c V^V 



• \^V 



Bl-Bl) 



^z\ 



\k x ky\ y \k x ky\ -\ (k y B y k x B x ) z 

e 2 2 

— ^B x B y z . 



(B21) 



Then, we classify all terms according to their order of 
iw/d< 1. 

The leading order is that of (w/d)^ 1 3> 1 and the off- 
block acquires the following main term 



H 



ih — * [k x (Jy + ky(T x ) k z 



(B22) 



It is important to remark that k x and k y contain only 
the z-component of the magnetic field. The transverse 
components B x and B y do not appear in H\h at this 
order of w/d. On this reason, the component B z has 
a larger effect on breaking the time-reversal symmetry 
than the other two components. 

The next order is that of (w/d)° ~ 1 and the off-block 
acquires the following correction 



H 



2m 



(*2 



.V373 



fc„ 2 ) + i {k x k y } a z 



-Vs^b [ «+ -.q 



ni 
2«73 



{zk z } ) (o x B x - a y B y ) . 
(B23) 



Further, there are two more orders: (w/d) 1 <?C 1 and 
(w/d) 2 <C 1, originating from terms containing z and z 2 , 
respectively. If they are required, one can find them by 
substituting Eq. (|B%t into Eq. ([Big]) . 
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